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Consider the following interacting particle system on the d -ary 
tree, known as the frog model : Initially, one particle is awake at the 
root and i.i.d. Poisson many particles are sleeping at every other 
vertex. Particles that are awake perform simple random walks, awak¬ 
ening any sleeping particles they encounter. We prove that there is 
a phase transition between transience and recurrence as the initial 
density of particles increases, and we give the order of the transition 
up to a logarithmic factor. 

1. Introduction. We study a system of branching random walks known 
as the frog model, and we discover a phase transition as the initial state 
becomes more saturated with particles. Similar phase transitions have been 
observed in related models, including activated random walk [9, 25], rein¬ 
forced random walk [17], killed branching random walk [1] and the contact 
process [18]. 

The frog model starts with a single particle awake at the root of a graph 
and sleeping particles at the other vertices. The initial configuration of sleep¬ 
ing particles can be deterministic or random. Particles that are awake per¬ 
form independent simple random walks in discrete time. When a vertex with 
sleeping particles is first visited, all of the particles at the site wake up and 
each begins its own walk. The name “frog model” was coined in 1996 by 
Rick Durrett; we continue the zoomorphism and refer to the particles as 
frogs. As with other interacting particle systems, the frog model is often 
motivated as a model for the spread of a rumor or infection (see [3], e.g.). It 
and its variants have also found interest as models of combustion [7, 22, 23], 
generally with particles moving in continuous time. 
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We call a realization of the frog model recurrent if the root is visited 
infinitely often by frogs and transient if not. Even if each individual frog is 
transient, the aggregate of visits to the root can still be infinite. For this 
reason, the transience or recurrence of the frog model gives a measurement 
of its growth, and the question of transience or recurrence for the frog model 
on a given graph is one of the most fundamental ones. 

The first ever published result on the frog model is that it is recurrent 
on T, d with one sleeping frog per site for all d [26]. In fact, the frog model 
on Z d is recurrent for any i.i.d. initial configuration of sleeping frogs [5]. It 
is natural to wonder if a sparser configuration changes the behavior. [20] 
exhibits a threshold at which a frog model with Bernoulli(a||a:||^ 2 ) frogs at 
each x € Z d switches from transience to recurrence. A similar phenomenon 
occurs when the walks have a bias in one direction: [14] finds that on Z, the 
model is recurrent if and only if the number of sleeping frogs per site has 
infinite logarithmic moment. Recently, this result was partially extended to 
Z rf in [10] and worked out in finer detail in [15]. 

Let Trf denote the full infinite d-ary tree, in which the root has degree d 
and all other vertices degree d+ 1. The question of transience or recurrence 
on is especially subtle. On one hand, the number of sleeping frogs grows 
exponentially with the distance from the root. On the other hand, each frog 
that wakes up has a drift away from the root; its probability of visiting the 
root shrinks exponentially as the starting vertex of the frog moves outward. 
The question of whether is transient for the one-per-site model is posed in 
[4] and again in [21] and [14], Surprisingly, the answer depends on the degree 
of the tree. In [16], we prove that the one-per-site frog model is recurrent on 
the binary tree and transient on d -ary trees with d > 5. 

We conjecture that the one-per-site frog model is recurrent for d = 3 and 
transient for d = 4. While we would like to pin this down and complete the 
picture of transience and recurrence for the one-per-site frog model on trees, 
we believe that the most interesting aspect of this work is that the frog 
model on trees is teetering on the edge between recurrence and transience. 
The point of this paper is to demonstrate this more precisely. We consider 
the frog model on with i.i.d. Poi(//) sleeping frogs at each site. Our result 
is a phase transition between recurrence and transience as p, varies. 

Theorem 1. Consider the frog model on a d-ary tree with Poi(/r) sleep¬ 
ing frogs per site. For all d> 2, there exists a critical value p c (d) > 0 such 
that the model is recurrent a.s. if p> p c (d) and transient a.s. if p < p c (d). 
The critical value satisfies 

Cd < pL c (d) < C'dlogd 
for some constants C and C'. 
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Proof. By a straightforward coupling, the probability of recurrence is 
monotone in //. By [16], Theorem 4, the probability of recurrence is either 
0 or 1. The theorem is then an immediate consequence of Propositions 6 
and 15, where we prove recurrence and transience, respectively. □ 

Contrast our result with the frog model on Z rf , which is recurrent for any 
i.i.d. configuration of sleeping frogs [5]. To show the existence of the recur¬ 
rence phase, we consider a restricted process that lets us take advantage of 
the recursive structure of T^. We then use a bootstrapping argument, show¬ 
ing that the number of returns to the root is stochastically larger and larger 
at each step. We establish the transience phase essentially by dominating 
the model with a branching random walk, using a similar argument as in 
[16]. As in that paper, the most difficult part is recurrence. Our result is 
an advance in that we are able to show recurrence on any d-ary tree with 
enough sleeping frogs. In [16], we prove recurrence only for d = 2, and the 
proof does not apply to a general choice of d; even extending it to d = 3 
seems difficult. The argument here relies on having Poisson many sleeping 
frogs at each site, however, and thus neither result implies the other. A more 
detailed comparison between the recurrence proofs in the two papers is in 
Section 2.4. 

Further questions. A nice general survey on the frog model can be found 
in [21]. Here, we pose four questions specifically related to the frog model 
on trees. 

The question most directly related to our paper is to better estimate the 
critical value /i c (d). We are interested in both the asymptotic behavior and 
precise values for small d. 

Open QUESTION 2. What is the correct order of /r c (d) as d-> oo? Also, 
what is the value of /i c (d) for small d? 

We suspect that /r c (d) = 0(d). As for the second question, the best bounds 
we can prove for d = 2 are 0.125 < fi c ( 2) < 1.13 (see Section 2.4). 

As a start at considering the frog model on less regular graphs, we would 
like to know if the analogue of our result holds on Galton-Watson trees. 

Open question 3. Consider a frog model with Poi(pi) frogs at each site 
of an infinite Galton-Watson tree. As /r varies, does a phase transition occur 
between transience and recurrence? 

We are also interested in the relationship between the frog model and the 
degree distribution of the tree. 
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Open QUESTION 4. Does the recurrence of the frog model on a Galton- 
Watson tree depend on the entire degree distribution or just the maximal 
degree? Concretely, consider a one-per-site frog model on a Galton-Watson 
tree where each vertex has probability p of having two children and probabil¬ 
ity 1 — p of having five children. [16], Theorem 1 implies that this is recurrent 
when p = 1 and transient when p = 0. Is it recurrent for any p <11 

This dependence on the maximal degree of the tree alone is seen in the 
contact process (see [18] and [19], Proposition 2.5). 

Our next question comes from Itai Benjamini and concerns the frog model 
on finite trees. Define the cover time to be the expected time for every frog 
to wake up in a one-per-site frog model on the full d -ary tree with height n. 
We call this the cover time since it is equivalent to the time when every 
site is visited. A naive bound on the cover time is 0{n 2 d n ), the expected 
time for a single random walk to visit every site, as shown in [2]. We have 
an unpublished proof improving this to 0(n 5 (d/V 2)”), but we suspect the 
correct value is polynomial. 

Open QUESTION 5. Is the cover time for the one-per-site frog model on 
a d -ary tree of height n polynomial in nl 

Possibly the cover time on finite trees relates to the recurrence and tran¬ 
sience properties on the corresponding infinite tree. For instance, it would 
be exciting to see that the cover time is polynomial in the height of the tree 
for d = 2 but exponential for higher d. This would be reminiscent of the 
contact process, which behaves similarly on finite lattices and trees as on 
their infinite counterparts [8, 11-13]. 

2. Recurrence. We start with a sketch. Let z/ be the law of the number of 
visits to the root in the frog model with Poi(/r) frogs at each site. To get some 
regularity, we restrict the motion of awakened frogs to the nonbacktracking 
component of their ranges. Call this the nonbacktracking frog model (more 
details are in Section 2.1) and let v be the law of the number of visits to 
the root in this model. A coupling argument in Proposition 7 confirms the 
intuition that 

(1) v<v. 

Here, A denotes stochastic dominance, that is u([x,oo)) <z/([x,oo)) for all 

x. 

In Section 2.2, we define an operator A under which the image of u has an 
interpretation in an even more restricted frog model. First, a bit of notation 
(see Figure 1) is necessary. Say the initial nonbacktracking frog moves down 
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Vd 



Fig. 1. The frog from 0 visits 0' and vi. Suppose at most one frog in the nonbacktracking 
frog model is allowed to enter each Td(fi) and only frogs woken at 0' and emerging from 
Td(vi) can enter other subtrees. We see in Lemma 9 that the number of visits to 0 is 
stochastically fewer than v and is distributed as An. 

the tree from the root 0 to 0 1 and then to v\. Let i> 2 .,Vd be the other 
children of 0' and let T d{vi) denote the subtree rooted at i>j. The measure 
Au is the law of the number of visits to the root in the nonbacktracking frog 
model with two further restrictions: 

(i) At most one frog can enter T d (vi) for each 1 < i < d. 

(ii) Only frogs woken at 0' and those emerging from TT^wi) can enter the 
other T d (vi). 

The advantage of (i) is that it makes the number of frogs emerging from the 
activated subtrees i.i.d. random variables. The advantage of (ii) is that it 
simplifies which subtrees become activated (see Lemma 8). Intuitively, these 
restrictions reduce the number of visits to the root. This is made rigorous 
in Lemma 9 where we prove that 

(2) Av < v. 

We stress that this a special property of v. In fact, the essence of our argu¬ 
ment is to show that when /r is large enough, (2) can hold only if v = 5^. 

Section 2.3 explores properties of A. In Lemma 10, we show that A is 
monotonic, meaning that for two probability measures it\ and 7T2, 

(3) if t\\ < 7T2, then Atti ■< AiT 2 - 

Lemma 11 shows that A acts nicely on the Poisson distribution. In fact, by 
writing the Poisson distribution in a nonstandard way (see Lemma 13), we 
can compare APoi(A) with Poi(A-fe). We carry this out in Proposition 14, 
where we show that when /r > 2 (d + 1) logd, there exists e such that 

(4) Poi(A + e) di A Poi(A) 

for all A > 0. This is where the value of fL plays a role. Proving (4) reduces 
to comparing two binomial distributions with parameters depending on fL. 

Now we explain how (1), (2), (3) and (4) imply the recurrence part of 
Theorem 1. 
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Proposition 6. If /a > 2(d+l)log<i, then the frog model is recurrent a.s. 
on the d-ary tree with an initial configuration of Poi(//) sleeping frogs per 
vertex. 

Proof. By (1), it suffices to prove that v is a point mass at infinity. 
From (2), we have 

Poi(O) A Av Y v. 

Statement (3) implies this relation is preserved under iterations of A. More¬ 
over, (4) lets us increase the Poisson term by e with each iteration. In sym¬ 
bols, this says that for all n > 1, 

Poi(en) Y A n v < A n ~ l v Y---YAvYv. 

Taking n —>• oo implies that v is a point mass at infinity, and so the frog 
model is recurrent almost surely. □ 

In the rest of this section, we will carry out this plan and prove state¬ 
ments (l)-(4). First, we give some notation. Recall that Y denotes stochas¬ 
tic domination. We also use the notation X Y Y to indicate that the law 
of X is stochastically dominated by the law of Y. An equivalent condition 
to stochastic dominance is that tt\ Y tt 2 if and only if there exists a cou¬ 
pling (X,Y) with X ~ 7Ti, Y ~ 7T2, and X < Y a.s. A thorough reference on 
stochastic domination is [24]. 

For a nonnegative random variable N, we use Poi(A^) to denote a mixture 
of Poisson distributions; when we write X ~ Poi(AI), we mean that X is 
coupled with N such that the distribution of X conditional on N = n is 
Poi(n). If N ~ 7r, we also use Poi(7r) to denote the same Poisson mixture. 
We similarly use the notation Bin(A r , jj) and Bin(7r,p). 

2.1. The nonbacktracking frog model. A random nonbacktracking walk 
on starting at a vertex xq moves in its first step to a uniformly random 
neighbor of xq. In all subsequent steps, it moves to a vertex chosen uniformly 
from all its neighbors except for the one it just arrived from. 

Suppose that ( S n ,n >0) is a random nonbacktracking walk starting from 
xq, stopped if it arrives at the root at step 1 or beyond. (If xq is the root, then 
it is never stopped.) Define the nonbacktracking frog model just as the usual 
frog model, except that the motion of a frog waking at xq is an independent 
copy of (S n ), rather than a simple random walk. The advantage is that when 
a nonbacktracking frog moves away from the root, it will forever remain in 
the just-entered subtree. This gives the model more self-similarity. As shown 
in [16], Proposition 7, (S n ) can be coupled with a simple random walk on 
I'd starting from xo so that its path is a subset of the simple random walk’s 
path. This lets us relate the nonbacktracking and usual frog models, proving 
(!)• 
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Proposition 7. Let v and z/ be the laws of the number of returns to 
the root in the nonbacktracking and usual frog models on T^, respectively, 
both with Poi (//) sleeping frogs per vertex. Then v -< v'. 

Proof. It suffices to show that we can couple the two models so that at 
least as many frogs visit the root in the usual model as in the nonbacktrack¬ 
ing model. We construct the coupling as follows. For each vertex v € T^. 
make the number of sleeping frogs on v identical in the two models. Make 
each frog’s path in the nonbacktracking model a subset of the correspond¬ 
ing frog’s path in the usual model as previously described. Thus, any frog 
woken in the nonbacktracking model is also woken in the usual model, and 
any visit to the root in the nonbacktracking model corresponds to a visit in 
the usual model. □ 

For the remainder of this section, we only consider the nonbacktracking 
frog model. We record an observation: Suppose the initial frog in the non¬ 
backtracking model steps from the root 0 to a child 0'. Since frogs are 
stopped at the root, no other child of the root besides 0' is ever visited, and 
all action occurs in the subtree rooted at 0'. 

2.2. Formal definition of A. Fix a probability measure n on the non¬ 
negative integers. We will define An to be the probability measure for the 
number of particles ending at 0 (see Figure 2) in the random system of 
nonbacktracking particles described below. 

The setting for the particle system is a star graph, consisting of a central 
vertex connected to d+ 1 leaf vertices. In a slight abuse of notation, we reuse 
the vertex names from Figure 1, calling the central vertex 0' and the leaves 
0 and ,Vd- Let X ~ Poi(/r) and X\,... ,Xd ~ n, all independent. Place 

X particles at 0' and Xj particles at each u,. Each particle if activated will 
perform an independent random nonbacktracking walk until it halts at a 
leaf. 

Initially, only the particles at 0' and at v\ are active. If one of these first- 
wave particles lands at Vi for i > 2, then the particles there are activated 
and begin independent nonbacktracking random walks until reaching a leaf. 
These second-wave particles do not activate other particles; only the first- 
wave particles have that power. The number of particles that finish at 0 in 
this system is a random variable, and we define An as its law. With these 
dynamics, we can summarize the system as follows: 

• Particles at 0' move to one of {0, v\ ,..., Vd} each with probability 1 /{d + 

!)• 

• Particles at v\ move to one of {0,V2, ■ ■ ■ ,Vd} each with probability 1/d. 

• If a first-wave particle visits Vi, the particles at V{ move to 0 with proba¬ 
bility 1/d. 


C. HOFFMAN, T. JOHNSON AND M. JUNGE 



Fig. 2. An interacting particle system, related to the frog model. Initially, the number 
of active particles at 0' is distributed as Poi(/x), and the number of active particles at 
v\ is distributed according to some probability measure n. Active particles take random 
nonbacktracking steps until reaching a leaf. For each 2 < i < d, if any of these particles 
reach n, then a new n-distributed batch of particles is released at Vi. These second-wave 
particles do not activate other vertices, (a) Initial state: particles at 0' and vi will move 
first and possibly release a second wave of particles from V2,...,Vd- (b) Terminal state: 
ff {particles at 0} ~ An. 


For 2 < i < d, let Ei be the event that a first-wave particle ends at v t . 
The following lemma follows from the definition of A. Informally, it says 
that conditional on how many of the events £ 2 , Eg occur, the number of 
second-wave particles ending at 0 is a sum of independent thinned copies 
of 7 r. 

Lemma 8. Conditional on 1 E t =u , the number of second-wave 

particles ending at 0 is distributed as the sum of u independent Bin(7r, 1 /d)- 
distributed random variables. 

Proof. If Ej occurs then by definition a 7r-distributed batch of particles 
is released at n*. With probability 1/d each released particle halts at 0. As 
particles move independently, the total number is distributed as Bin(7r, 1/d). 
Since the second-wave particles cannot wake other sites, the total number 
of particles to arrive is distributed as claimed. □ 

Now, we show the connection between this operator and the frog model. 

Lemma 9. Let v be the distribution of number of returns to the root in 
the nonbacktracking frog model on the d-ary tree with sleeping frog distribu¬ 
tion Poi(/r). Then Av -<v. 
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Proof. Let T t /(x) denote the subtree of rooted at a given vertex x. 
Recall that no children of the root other than 0the child visited by the 
initial frog, are ever visited. In light of this, it will be helpful to think of the 
nonbacktracking frog model as taking place on 0 U T d (0') rather than on 
all of T d . 

We say that the frogs sleeping on some vertex v G Td(ui) wake within 
Trf(ui) if there exists a chain of vertices xi,...,x m = v all in T^ui) such 
that the initial frog starting from the root visits xi, a frog starting at x\ 
visits X 2 , and so on. More simply, a frog is woken within Trf(ui) if it would 
have been woken even if there were no frogs sleeping on any vertices outside 
of T d (vi). 

We define some random variables counting frogs that might possibly visit 
the root. Let X ~ Poi(/x) be the number of frogs sleeping on 0', which are 
woken by the initial frog. Let X\ be the number of frogs waking within 
T^(ui) that visit 0'. We claim that X\ is distributed as v. Indeed, when 
we consider frogs as waking only if they wake within Trf(ui) and relabel the 
vertices {0'} U T^(ui) as {0} U T d (0'), we see a process identical in law to 
the original nonbacktracking frog model. Call the frogs counted by X and 
X\ the first-wave frogs. 

For each 2 < i < d, let E t be the event that some of the frogs counted 
by X or X\ move to u,. Conditional on E^. arbitrarily choose one of these 
frogs that visits Vi and call it /. We say that the frogs at v are woken 
within Trf(uj) if there exists a chain of vertices xi,...,x m = v in T d (vi) 
such that / visits x \, a frog starting at x\ visits X 2 , and so on. Let X % 
be the number of frogs waking within T^(uj) that visit 0'. By the same 
argument showing that X\ ~ n, the distribution of X t conditional on Ei is 
also v. Furthermore, for any {R,..., i C {2,..., d}, the random variables 
X n ,..., X lk are conditionally independent given E ^,..., E lk , since each X t 
is determined solely by the paths of the frogs sleeping in T d (vi). We call the 
frogs counted by X 2 , ■ ■ ■ ,X d the second-wave frogs. 

The first- and second-wave frogs all visit 0 r . We define V" as the number 
of these that move from there to 0. 

Claim. V"^An. 

Proof. Our strategy is to show that the first-wave frogs behave iden¬ 
tically as the first-wave particles, and then to show that the second-wave 
frogs conditional on the behavior of first-wave frogs behave the same as the 
second-wave particles conditional on the behavior of the first-wave particles. 

For the first of these claims, consider the first-wave frogs, counted by 
X and X]. Observe that X and X\ are independent with X ~ Poi(/r) and 
X\ ~ v, just as in the particle system defining An. The frogs counted by X 
move from 0' independently to a random choice out of 0,v±,v d , and the 
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frogs counted by X\ move from 0' independently to a random choice out of 
0, V 2 , ■ ■ ■, Vd, also matching the particle system. Thus, the locations of the 
first-wave frogs one step after leaving 0' are distributed identically to the 
ending locations of the first-wave particles. 

Now, condition on some arrangement of the first-wave frogs on 0, v \,..., Vd 
one step after leaving 0. Suppose that u out of the vertices V 2 ,,Vd are oc¬ 
cupied by first-wave frogs in this arrangement. The number of second-wave 
frogs visiting 0' conditional on this arrangement of first-wave frogs is a sum 
of u independent copies of v. Each second-wave frog that visits 0' has an 
independent 1/d chance of moving next to 0. Thus, the number of second- 
wave frogs that visit 0 is the sum of u independent copies of Bin(z/, 1 /d). 
This matches the conditional distribution of second-wave particles ending 
at 0 given in Lemma 8. Thus, the distribution of the number of first- and 
second-wave frogs visiting 0 is the same as the distribution of the number 
of first- and second-wave particles ending at 0, which is by definition Au. 
□ 


With this claim, the proof of the lemma is almost complete: Let V be the 
total number of visits to 0 in the nonbacktracking frog model. Since V" < V 
with V" ~ Av and V ~ v, we have shown that Av -< v. □ 

2.3. Properties of A. We first show (3), monotonicity of A with respect 
to stochastic dominance. 

Lemma 10. If ir i ^ 7r2, then An\ P An 2 ■ 

Proof. If tt\ < 7r2, then we can couple the two particle systems defin¬ 
ing Ajt\ and Ait 2 so that the second particle system contains all the same 
particles as the first, moving identically, as well as additional ones. Thus, 
at least as many particles visit 0 in the second system as in the first, and 
An 1 -< Att 2 - □ 

Now, we describe the result of applying A to a Poisson distribution, whose 
thinning property simplifies things. 

Lemma 11. The distribution ^lPoi(A) is a mixture of Poisson distribu¬ 
tions, given by 

(5) ^P 0l( A H P„i((k±ih + _^ T ), 




where 

( 6 ) 


U ~ Bin d — 1,1 — exp 
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Proof. In the particle process defining AlPoi(A), let Y u ^ v be the num¬ 
ber of particles that start at u and finish at v, for u£ {0',vi,...,Vd} 
and v € {0,vi,... ,Vd}. Each of the Poi (p) particles starting at 0' moves 
to a random neighbor. By Poisson thinning, the random variables Y 0 /_ > . v 
for v £ {0,vi, ... ,Vd} are independent and distributed as Poi(^/(d + 1 )). 
Similarly, Y V1 _ >v for v € {0, V 2 , ■ ■ ■, Vd} are independent and distributed as 
Poi(A/d). These two collections of random variables are also independent of 
each other. 

Thus, the number of first-wave particles that move to u* for each 2 < 
i < d are independent and distributed as Poi(A/d + p/(d + 1)). Let U be 
the number of vertices out of {v 2 ,...,Vd} that are visited. As each vertex 
has an independent 1 — exp(— X/d — p/(d + 1 )) chance of being visited, the 
distribution of U is as given in (6). And since U is determined by Y 0 t^ v . 
and Y Vl - ¥Vi for i = 2, ..., d, it is independent of Y 0 i^ 0 and Y Vl ^ 0 . 

By Lemma 8 and Poisson thinning, the number of second-wave particles 
ending at 0 is Poi(t/A/d). The number of first-wave particles ending at 0 is 
Y 0 r_f 0 + Y V1 _> 0 , independent of U and distributed as Poi(A/d + fi/(d+ 1 )). 
Summing these together yields (5). □ 

We are nearly in a position to establish that A n Poi(O) grows without limit 
as n —>• oo. First, we need two technical lemmas on the Poisson distribution. 


Lemma 12. Let Z\ be distributed as Poi(A) conditioned to be nonzero. 
If Ai < A 2 , then Z Xl Y Z X2 . 


Proof. Consider the Radon-Nikodym derivative of the law of Z X2 with 
respect to the law of Z Xl , 



The function r(k) is increasing, and it is straightforward to show that this 
implies that Z Xl Y Z X2 (or see [24], Theorem l.C.l). □ 


Lemma 13. Let Z ( x / n , Z X j n ,... be independent and distributed as Poi(A/n) 
conditioned to be nonzero. Let M be independent of these and be distributed 


7(1) 7(2) 

^A/n’^A/n’ ' • 


as Bin(n, 1 — e x ^ n ), and let 


M 



i— 1 


Then Z is distributed as Poi(A). 
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Proof. Decompose Poi(A) as a sum of n independent copies of Poi(A/ro). 
Let M be the number of these that are nonzero, and condition on M to get 
the desired representation. □ 

Finally, we prove (4). 

Proposition 14. If /r > 2 (d + 1) logd, then there exists s > 0 such that 

Poi(A + e) < AlPoi(A) 

for all A > 0. 


Proof. Let X ~ Poi(A + e) for some e > 0 to be chosen later, and let 
Y ~AlPoi(A). We start by decomposing X into a sum of Poissons condi¬ 
tioned to be nonzero. For any a, let be distributed as Poi(a) 

conditioned to be nonzero, and let Z a ~ Poi(a) (with no conditioning). Take 
all these random variables to be independent. By Lemma 13, we can write 
X as 

M 

m * = z (x+c)/i +Y,ztL)M< 

1=1 


where 


M ~ Bin 


d — 1,1 — 


exp 


A T £ 
d 


( 8 ) 


We now turn to Y, which by Lemma 11 is distributed as 

(U + 1)A 


Poi 


d 


+ m 


where m = fj,/(d + 1) and U ~ Bin(d — 1,1 — exp(— X/d — m)). Let Y' ~ 
Poi((f7 + 1)(A + m)/d). For each u, the distribution of Y' conditional on 
U = u is stochastically dominated by the distribution of Y conditional on 
U = u, simply because Poi(a) ^ Poi(6) when a <b. It follows that Y' -< Y. 
Thus, it suffices to show that X ^ Y'. Decomposing Y' by Lemma 13 and 
using the same notation as before, we can write Y' as 


( 9 ) 


N 

* ' = Z (\+m)/d + ^2 Z {\+ m )/d 
2=1 


with 


N ~ Bin 


(u, 1 - 




exp 


A + m\ A 

~ir JJ 
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These decompositions allow us to stochastically compare X and Y'. As¬ 
sume that e is chosen to be smaller than m. We claim that to show that 
X ^ Y', it suffices to show that M < N. Indeed, we can then couple the 
random variables on the right-hand sides of (7) and (9) so that: 


(1) M < N] 

( 2 ) Z(\+e)/d — Z(\+m)/di 

(3) Z^ x+£)/d < Z$ +m)/d for each i. 

Property (2) is possible because Poi(a) ^ Poi(6) if a<b, and (3) is possible 
by Lemma 12. Together, this yields a coupling of X and Y' with X < Y'. 

Thus, it only remains to show that M -< N. Recalling that U is itself 
binomial, we have 

N ~ Bin ^Bin — 1, 1 — exp ~ ^ > 1 — ex P ~d~ 



= Bin 




1 — exp ( — — — m 


1 — exp — 


A + 771, 

~d~ 


Since M and N are both binomial, proving M < N reduces to comparing 
their parameters. The argument will be complete once we show for some 
e > 0 and all A > 0, 


do) l-exp(-MT) 


< 


^1 — exp 



1 — exp 


A + m\ A 

)) 


Some basic calculus (see Lemma 16 in the Appendix) establishes that for 
all d> 2, 


g—2 log d _|_ g—2 log d/ d < 


Since m > 2 log d, we can choose e > 0 such that 

l>exp^-0 > e - m + e- m / d . 
Multiplying both sides of this inequality by e~ x ^ d gives 


exp 


A T £ 


> exp ( — — — m ) + exp ( — 


V d 


A + m 


\ d 


Thus, 


i ( A + e Wi ( A \ l 

1 — exp--— < 1 — exp —- — ml— exp — 


V d 


A + m 


V d 


A 


< ( 1 — exp ( — — — m ) ) ( 1 — exp ( — 


A + m 


Looking back at (10), we have shown that M < N. □ 


We have now proven (l)-(4), completing the proof of Proposition 6. 
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2.4. Comparison to one-per-site results. In [16], we proved that the frog 
model on a binary tree with one sleeping frog per site is recurrent. The proof 
has the same overarching idea as here: We use the self-similarity of the tree 
to obtain a recursive distributional relationship for the number of returns 
to the root. We then use this relationship in a bootstrapping argument, 
assuming that the number of visits to the root is stochastically larger than 
Poi(A) and proving that it is in fact stochastically larger than Poi(A + e). 

The major difference between the two arguments is in the bootstrapping 
portion. The approach in this paper using traditional stochastic domination 
fails with the one-per-site frog model. The problem is that the distributions 
given by successively applying the analogue of the A operator in the one- 
per-site model have finite support, and hence are never stochastically greater 
than any Poisson distribution. Our proof in [16] instead uses an exotic defini¬ 
tion of stochastic dominance, where 7Ti is dominated by Ti 2 if the probability 
generating function of tt\ is greater than the probability generating function 
of 7r 2 . 

This generating function approach works better than the technique in 
this paper in some ways and worse in others. On one hand, it can handle 
both deterministic and random initial configurations. On the other hand, 
the generating function approach seems confined to small values of d. It 
relies on a purely analytic argument that is elementary but difficult. It seems 
impossible to apply this argument to an arbitrary choice of d. Even for d = 3, 
the generating functions to be analyzed become extremely complicated. The 
technical advance in this paper is the probabilistic argument we give in 
Proposition 14, which allows us to work on any d- ary tree. 

3. Transience. The main idea of our proof of transience is to consider 
a weight function on the frog model. To analyze the weight function, we 
bound the frog model by a branching random walk. The weight function is 
the frog model analogue to a common martingale derived from branching 
random walk (see [6]). 

Proposition 15. If Eg < y u , then the frog model with an indepen¬ 
dent copy of rj frogs per site on is almost surely transient. 

Proof. Let F n be the set of frogs awake at time n. For / € F n , let |/| 
denote the level of / on the tree (that is, its distance from the root). We 
define a weight function 



f&Fn 

with 6 to be chosen shortly. Let 
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Before we explain the meaning of this, we minimize m by setting 9 = 
log((E/y + l)d)/2, making 

2\/ (E?? + 1 )d 

m = — 1 - —-— -< 1 

d+1 

under our assumption that E?/ < ^ ^ ■ 

The strategy of the proof now is to show that W n — > 0, and hence that the 
root eventually stops being visited. The term m gives an upper bound for 
the expected contribution to W n +1 of a frog at time n in the following way: 
Suppose that at time n, some frog / is at level i of the tree for any i > 1. 
With probability l/(d+ 1), the next jump of / is toward the root, waking no 
frogs. With probability d/(d + 1), the jump is away from the root, possibly 
waking up an ry-distributed number of frogs. Thus, the expected contribution 
to W n +1 from / and any frogs it wakes at time n + 1 is at most e~ Sl m. If 
/ is at the root at time n, then the expected contribution to W n +1 from / 
and the frogs it wakes is at most E[?y + l]e~ e , which is bounded by m given 
our choice of 9. Therefore, 

E[W n+ i \W n ]<Y, e" 0|/l m = mW n . 

f€F u 

Thus, W n /m n is a positive supermartingale. By the martingale convergence 
theorem, it converges almost surely to a finite limit. Since m n —>• 0, we also 
have W n —>• 0 a.s., which implies that eventually no frogs are present at the 
root. □ 


APPENDIX 

Lemma 16. x~ 2 + x~ 2 ^ x < 1 for all x>2. 

Proof. Let f(x) = x~ 2 + x~ 2 / x . First, we show the inequality holds on 
the interval [2,8]. Since x~ 2 is decreasing, 

f(x)< \ + x~ 2/x . 

It is easily checked that the maximum of x~ 2 ! x on [2,8] occurs at x = 8 and 
is less than |. 

Next, we consider x>8. L’Hopital’s rule implies that lirn^oo/(x) = 1. 
Thus, it suffices to confirm that f{x) is increasing on [8, oo). We compute 

f(x) = 2x- (2 / x) - 2 (logx - x^- 1 - 1). 

For x > 8, it holds that x ( 2 /*) -1 < 1. Hence, 

f{x) > 2 x~ ( k2 / x ' , ~ 2 (\ogx — 2 ), 

which is positive on [8, oo) since log8 >2. □ 
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